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Abstract 

We use the leading singularity technique to determine the planar three-loop five-particle ampli- 
tude in N = 4 super Yang-Mills in terms of a simple basis of integrals. We analytically compute 
the integral coefficients for both the parity-even and the parity-odd parts of the amplitude. The 
parity-even part involves only dual conformally invariant integrals. Using the method of obstruc- 
tions we numerically evaluate two previously unfixed coefficients which appear in the three-loop 
BDS ansatz. 
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I. INTRODUCTION 



Scattering amplitudes in gauge theories are remarkable objects with many properties 
hidden in the complexity of their Feynman diagram expansion. It is natural to expect the 
most symmetric gauge theory, Af = 4 super- Yang-Mills (SYM), to be especially rich. Two 



very recent papers [1 



2j have used AdS/CFT to relate and shed light on two particularly 
remarkable properties of scattering amplitudes in SYM. These properties are dual confor- 
mal invariance, first observed for MHV amplitudes in [3] and conjectured to extend to all 
amplitudes in J4], and an equality between amplitudes and certain lightlike Wilson loops, 



first established at strong coupling in and at weak coupling in 



were pushed Jorward in a number of papers including 
reviewed in 



10 



11 



3L 



12 



13 



hese developments 
"1 

14] and have been 



15} , where a comprehensive set of references may be found. 



Much of the recent interest in multi-loop scattering amplitudes has been stimulated by 



the ABDK/BDS ansatz 16|, ll7|] which suggested that multi-loop MHV amplitudes satisfy a 



powerful iteration relation implying a simple exponential form for the full all-loop amplitude. 
Although the ABDK/BDS ansatz was successfully tested for four particles at two 16( and 



three 



21 



17l | loops, as well as for five particles at two loops 



18 



191 ] . some doubts raised in [20, 



221 ] necessitated an explicit calculation of the two- loop six-particle amplitude 11] which 
conclusively demonstrated the incompleteness of the BDS ansatz. 

Indeed six particles is the earliest that the hypothesized dual conformal symmetry of 
amplitudes could have allowed BDS to break down: for n = 4, 5 the symmetry fixes the form 



2p| . Beginning at n 



6 dual conformal 



of the amplitude up to a few numerical constants [9 ! , 
invariance becomes substantially weaker, determining the form of the amplitude only up to 
an arbitrary function of dual conformally invariant cross-ratios. Nevertheless it was found 
in parallel work 



11 



121 ] that the amplitude/Wilson loop equality survives at two loops for 



n = 6 despite not being required by dual conformal invariance. 
In this paper we study the three-loop BDS ansatz 



{e) + \{M^{e)f-M^\e)M^( 



/ (3) (e)M«(3e) = C (3) + 0(e) 



where is a previously undetermined numerical constant. In sections II and III we use the 

n 

leading singularity method [23] to determine the (four-dimensional cut-constructible part of 
the) 3-loop 5-particle amplitude in terms of a simple basis of integrals. In section IV we 
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then numerically evaluate enough pieces of these amplitudes (the pieces called 'obstructions' 



in 



24 



251 ] ) to determine = 17.8241. Although current developments strongly suggest 



that the quantity appearing on the right-hand side of (II. ip will in general be non-constant 
(but still dual conformally invariant) for n > 5, there is some utility in knowing the precise 
number = 17.8241 since for any n, whatever appears on the right-hand side of (II. ip 
must approach this same number in any collinear limit. 



II. OUTLINE OF THE CALCULATION 



Our goal is to find a compact expression for the planar 3-loop 5-particle amplitude in 
Af = 4 SYM as a linear combination of some basic integrals. Several powerful and related 



techniques 



27 



28 



29. 



or carr ying out calculations such as these include unitarity based methods 



ad, 



31 



32J and more recently, building on 33j, |34j], maximal cuts 35[] and the 



leading singularity method [23[. For the present calculation we find it convenient to use the 
leading singularity method (see also [36|, |37j) since it allows for all integral coefficients to 
be determined analytically by solving simple linear equations. In this section we provide a 
detailed outline of the steps involved in setting up the calculation. 



A. Review of the Leading Singularity Method 

Suppose we are interested in calculating some L-loop scattering amplitude A. On the 
one hand, the amplitude may of course be represented as a sum over Feynman diagrams Fj, 

A(k)=J2 I WdXF^l), (2.1) 

j J a=l 

where k are external momenta and £ a are the loop momenta. However it is frequently the 
case, especially in theories as rich as Af = 4 SYM, that directly calculating the sum over 
Feynman diagrams would be impractical. Rather the calculation proceeds by expressing A 
as a linear combination of relatively simple integrals in some appropriate basis {/, }, 

A(k) =J2°i( k ) I H dd£ - J *(^)' ( 2 - 2 ) 

i J a=l 

and then determining the coefficients q by other means. 
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With the leading singularity method we equate (12. II) and (12. 2p and perform the integral 

Vc;(A;) f d 4 £Ii(k,£)= [ AVFjfM) (2.3) 

over contours V G C 4L other than the real £-axis. At L loops each contour is a T AL inside 
C 4L . For each contour T we obtain one linear equation on the coefficients q. Of course 
if T is a random contour then we would generally get the useless equation = 0, so we 
should choose contours such that the integral on the right-hand side of (I2.3P evaluates the 
residue on the isolated singularities of Feynman diagrams, which are associated with the 
locus where internal propagators become on-shell. Since the number of isolated singularities 
in a generic L-loop diagram can be as high as 2 L (simple diagrams can have fewer isolated 
singularities), the leading singularity method gives rise to an exponentially large (in L) 
number of linear equations for the coefficients q. We note that the homogeneous part of 
these linear equations (the left-hand side of (12. 3p ) depends only on the set of integrals {/j} 
and the choice of contours, while the details of which particular helicity configuration is 
under consideration enters only into the inhomogeneous terms on the right-hand side. 



B. Integration Strategy: Collapse and Expand 



Here we briefly review from 23|, |34], |35|, |36|, |37| the integration rules which make it simple 



to evaluate the contour integrals appearing in (12.31) in the cases relevant to the present 
calculation. Let us focus on a box with loop momentum p and external momenta fcj. The 
box may be sitting inside a higher-loop diagram, in which case the ki may involve other loop 
momenta. The sum over Feynman diagrams contains poles at the locus 

S = { P EC 4 :p 2 = 0, (p-Aa) 2 = 0, (p-k 12 ) 2 = 0,(p + h) 2 = 0}, (2.4) 

which, for generic ki, consists of two distinct points. To each of these points there is an 
associated contour T p such that integrating p over T p calculates the residue at the associated 
point. 

The residue of a one-loop amplitude at one of these poles is computed by removing the 
four internal propagators and evaluating the product of on-shell tree amplitudes at the four 
corners (summed over all helicities of internal states). In the simplest application, when all 
four ki satisfy kf = 0, this product evaluates on either contour T p to a four-particle tree 



amplitude, leading to the 'collapse rule' graphically depicted as 




The figure on the left indicates the sum over that subset of all one-loop Feynman diagrams 
in which all four of the indicated propagators are present. Of course it may as well be the 
sum over all one-loop Feynman diagrams since those that do not contain all four of the 
indicated propagators contribute zero to the residue. 

When one of the kf is non-zero, the result (12.51) holds on only one of the two T p contours, 
while the integral over the other contour gives zero. Given a helicity assignment for the 
external particles it is a simple matter to determine which of the two solutions leads to the 
non-zero result. 

It is frequently the case that after collapsing a box in some loop momentum p there are 
less than four exposed propagators in some other loop momentum, which would apparently 
indicate a codimension 1 singularity rather than an isolated singularity. In this case one can 
use the 'expand rule' 




+ terms non- singular at {ki + k 2 ) 2 = 



k 2 h 



V 

= v + terms non-singular at ik 2 + k 3 ) 2 = 

A 

ki fc 4 

to expose additional propagators inside a tree amplitude. The choice of how to expand 
is correlated with the choice of integration contour for the next loop momentum. In the 
example shown here, the terms isolated on the first line are those which survive a contour 
integration around the singularity at (ki + k 2 ) 2 = while the second expansion displays 
those terms isolated by a contour integration around the singularity at (k 2 + k^) 2 = 0. 



5 



These two simple rules are sufficient for evaluating all contour integrals appearing on the 
right-hand side of (12. 3p in this paper. Finally, scalar integrals appearing on the left-hand 
side are integrated via the simple rule 

4 1 1 

r p dp p 2 (p- h) 2 (p - h 2 ) 2 (p + h) 2 = (h + k2) 2 (k 2 + ' (2 ' 6) 

which is valid as long as at least three of the hi satisfy kf = (we will not encounter any 
other cases in the present calculation). We have chosen a simple normalization factor of 
1 on the right-hand side of (12. 6ft ; this will be adjusted below in (13.31) to match standard 
conventions for normalizing amplitudes. 

C. Choosing a Sufficient Set of Contours 

In order to proceed systematically we begin by enumerating all planar 3-loop 5-particle 
topologies which are free of tadpoles, bubbles, and triangles, since such diagrams are un- 



neccesary due to TV = 4 supersymmetry (see 39( for a thorough discussion). This leaves 17 



topologies, of which 5 do not have any associated leading singularities and are therefore of 
no interest to us. The remaining 12 topologies are shown in fig. [IJ 

Each topology in fig. [1] has several distinct associated leading singularities, each of which 
gives rise to an equation via (12.31) . The information contained in this collection of equations 
is highly redundant — the equations obtained from only a small subset of the leading sin- 
gularities are sufficient to determine all coefficients, while the remaining equations serve as 
consistency checks. We now present a few details explaining how to extract a set of equa- 
tions sufficient for determining all coefficients. We have verified a number of the additional 
equations to check consistency, but have not performed an exhaustive search for all possible 
leading singularities. 

The topologies fall naturally into three different categories according to how we choose to 
implement the collapse and expand rules. Let us now address each category in turn, giving 
in each case the details of the simplest topology as an example. 
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FIG. 1: The planar 3-loop 5-particle topologies associated to leading singularities. Each figure 
represents a sum over that subset of Feynman diagrams in which all of the indicated propagators 
are present. We label the external momenta clockwise with k\ at the leg indicated with the arrow. 
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1. Example 1: Topology L 



Topology L has several leading singularities, but the simplest ones can be isolated as 
indicated in the following cartoon: 




In words: we first integrate the sum of Feynman diagrams over a p contour which collapses 
the associated massless box, then expand around (q + ki) 2 = keeping only the singular 
terms indicated. Integrating q over an appropriate contour isolates these singular terms 
while collapsing the massless box. The final integral over r is again accomplished using the 
collapse rule. 

The leading singularities exposed by these steps are those located at the locus 

S L = {(p,q,r) EC 12 : p 2 = 0, q 2 = 0, r 2 = 0, (p + k x f = 0, (r - k 5 ) 2 = 0, 

(r - k i5 ) 2 = 0, (r + k l2 f = 0, (q+ k l2 f = 0, (p - k 2 ) 2 = 0, 
(q - r) 2 = 0, (p + q + h) 2 = 0, (q + h) 2 = 0}. (2.8) 

For generic external momenta the set Sl consists of 8 distinct points in C 12 . For each point 
in Sl there is an associated contour which computes the residue at the point and hence leads 
to an equation via (12. 3p . 

It remains only to construct an appropriate ansatz for the left-hand side of ( 12.31) . We try 
a linear combination of the two most natural integrals of topology L, 




Here and in what follows we use pictures as shorthand for the corresponding scalar inte- 
grands, so for example the first term in (12. 9p represents 



^ p 2q2 r 2(p _|_ fc 1 )2( r _ k 5 ) 2 (r - k i5 ) 2 (r + k 12 ) 2 (q + h 2 ) 2 (p - k 2 ) 2 (q - r) 2 (p + q + h) 2 

(2.10) 
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while the dotted line in the second picture in (12.91) indicates a factor of (r + ki) 2 in the 

numerator of the integrand. 

Integrating (12. 9p over the contours detailed above leads to the expression 

L + L^r + h) 2 
sj 2 s u s 4 5(r + fci) 2 ' 

where the denominator factors arise from the Jacobians in (12.61) . Equating this to the result 
of (12.71) and choosing a particular helicity configuration leads to the equation 

itY r (r + +f)* = Aree(1 ~' 2 ~ ' 3+ ' 4+ ' 5+) ^ (2 ' 12) 

S 12 S3 4 S45(r + K\) 

Of course this must be evaluated on the locus Sl, and it is easy to check that Sl contains 
only two different values of r: 

n = A 5 ( A 5 + ^4vA 4 ^) , r 2 = ( A 5 + A 5 , (2.13) 



(5,3) V V [5,3]' 

giving us two distinct equations which are sufficient to determine the coefficients L and L\ 
uniquely. 

Topologies D, G, and iV proceed in exactly the same manner, except that in these cases 
more than two integrals appear on the left-hand side. Topologies A, C, E and K are 
also very similar, except that since these three topologies only have 10 exposed propagators 
(rather than 11) it is necessary to isolate a second hidden singularity by performing a second 
expansion prior to integrating over r. 



2. Example 2: Topology M 

For topology M it is sufficient to consider even simpler contours. We first collapse and 
expand the p box as done above for topology L, arriving at 



(2.14) 



r w q 

At this stage it is convenient to integrate over a symmetric contour of the type considered 




m 



231 ] where we require that (q, r) and [q, r] separately vanish instead of just (q + r) 2 = 0. 



This leads us to consider the locus 



S M = {{p, q, r) G C u : p 2 = 0, q 2 = 0, r 2 = 0, (p - h) 2 = 0, (r + h) 2 = 0, 



(q + h 5 ) 2 = 0, (r + k l2 f = 0, (p-q- h 5 ) 2 = 0, 
(p-k 45 ) 2 = 1 (g + A; 5 ) 2 = 0, (q,r) = 0, [g,r] = 0}. (2.15) 



For generic external momenta Sm consists of 4 isolated points, each of which leads to one 
linear equation for the integral coefficients. Note that the right-hand side of (I2.3P always 
vanishes for such symmetric contours since the associated product of tree amplitudes must 
vanish when (q, r) = = [q, r\. 

For topologies B, J and G we proceed along exactly the same lines (we already treated G 
in the first example, but additional equations are needed to fix all of the coefficients which 
appear for this topology). It turns out that for topology J an interesting and very useful 
feature emerges: here we model the left-hand side as the linear combination 



We will not display all of the relevant integrals explicitly, but they all have the property 
that they either vanish on the locus Sj (so that they do not enter the associated equations), 
or they contain the same numerator factor as J\ shown here, which we denote by £ 2 . Now 
when we perform the q integral one of the Jacobian factors is l/£ 2 , so we obtain the equation 



Since £ 2 = on the locus Sj, we immediately see that the coefficient J must vanish in order 
to avoid a contradiction. Perhaps a safer way to express this is to say that we can consider 
an equation obtained by multiplying both sides of (12.31) by £ 2 before performing the contour 
integrals. Having determined that J = 0, we then see that (I2.17P gives an equation relating 
J\ to the other coefficients. This trick is also useful for other topologies, in particular for C 
and E. 

3. Example 3: Topology R 

There are five different triple-box 9-propagator topologies, of which topology R is the 
only one with associated leading singularities. These are situated on the locus 



S R = {(p,q,r)eC 12 : p 2 = 0, (p + q + k 15 ) 2 = 0, q 2 = 0, (q - h) 2 = 0, (p - r) 2 = 0, 




(2.16) 



— + J\ + several other coefficients = 0. 



(2.17) 
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(r - k 23 ) 2 = 0, r 2 = 0, (r + q + k 15 ) 2 = 0, (p + A^) 2 = 0, 
(g + A; 15 ) 2 = 0, (r + A^) 2 = 0, (r + A; 15 ) 2 = 0}. (2.18) 

Here the first nine conditions are the visible propagators, while the last three are hidden 
singularities. In order to see what the right-hand side of (12.31) should be let us begin by 
integrating out q to collapse the first box. This leads to 

(2.19) 

For the first time we find a triangle-triangle diagram rather than a triangle-box or box-box. 
The Jacobian factor from integrating the corresponding scalar integral is l/(q+ki5) 2 (r + ki) 2 , 
suggesting that we expand (12.191) to expose either l/(q + &15) 2 or the l/(r + k\) 2 propagator, 
but it is clearly impossible to expand both simultaneously. Either choice leaves us with a 
sum over triangle-box Feynman diagrams, which vanishes due to M = 4 supersymmetry. 
The right-hand side of (I2.3P is therefore zero for the R topology leading singularities in 
eq. (12TTB1 . 




III. THE 3-LOOP 5-PARTICLE AMPLITUDE 



A basis of integrals which is sufficient for representing all of the leading singularities of 
the amplitude is shown in fig. [2j By solving the collection of linear equations as explained 
in the previous section we find their coefficients 

A 2 1-73 

73-73 

R 75 
-D — — S12S23S34S45- 



75 - 75 ' 

r - , , 2 73 



" 73 - 73 
Di = S12S23S51 — , 

73 - 73 

1-74 

L>2 = -S23S34S15 — , 

74 - 74 

& = — S12S23S51 — , 

73 - 73 

Gl = Si2S2 3 S 51 — , 

73 -J3 
72 

G*2a = — •S23S45S51 — , 

72 - 72 
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FIG. 2: The 17 independent integrals appearing in the ansatz. Other integrals can be obtained 
by rotations or reflections. As in fig. Q] we label the external momenta clockwise with k\ at the 
position indicated by the arrow. 



G2b 

Ji 
K 
L 

u 

M 
N 
N, 
R 

where we have introduced the quantity 

li = 



1 



— S23S34S45- 

75 - 75 

1 

— S34S45S51 — . 

7i - 7i 

3 1-73 

S 12 S23 — , 

73-73 
1 



73 ~J3 

2 71 

— S 12 S34S45 — . 

71 - 7i 
2 1 

— s 12S 45 S5i — . 

72 - 72 



-S51-Sl2S34S 45 - 



7i_- 7i 
7i 



<S23S45S51" 



71 - 71 

1 - 71 

7i - 7i 



(z + 2,i + 3)[i + 3, 
(i + 2,i + 4)[i + 4,i\ 

and 7 is given by the parity conjugate of this expression (i.e., (a, b) 



(3.1) 



(3.2) 



[a, b]). In each case we 



have suppressed an overall factor of AI e e . In order to connect to the standard normalization 



12 



conventions used in the study of the BDS ansatz it is necessary to multiply by an overall 
factor of (— 1/2) L . The complete amplitude is therefore assembled via the formula 



M, 



(3) 



A®/Af" = -- ^coefficient^ x integral;. 



(3.3) 



permutations integrals 

The first sum runs over the 10 cyclic and anti-cyclic orderings of the labels 1, 2, 3, 4, 5 of the 
external particles and the second sum runs over the 17 integrals in fig. [2j Si is a symmetry 
factor to compensate for possible overcounting: S = 2 for integrals B, Di, D 2 , L, Li, 
N, and R, and S = 1 for the others. 

The presentation (13.11) makes it simple to read off the parity-even parts of the coefficients, 
which will be useful in the following section. We find the parity-even parts 



A i = 2 S i2S 23 , 



G 



2d 



2 S 23S45S51, 



B = -^Si 2 S 2 3S34S45, 
D2 = 2 S 23 s 34"S5l, 
K = J s 12 s 23, 



L\ — -S12S34S45, 



-S12S34S45, 



R = --S23.s45.s51- 



(3.4) 



with all others vanishing. We note that only the coefficients associated to dual conformal in- 
tegrals have non- vanishing par ity-even parts, as expected based on the pattern of previously 



studied amplitudes [U 



17 



19 



39|. 



IV. THE THREE-LOOP BDS ANSATZ 



itudes in gauge theory are very 



The infrared divergences of higher loop scattering amp 
simply related to those of lower loop amplitudes |40j]. In 
in M = 4 SYM this simplicity persists, at least for MHV amplitudes, to the finite terms 



17j |. it was conjectured that 



as well. Although the ex 
earlier doubts raised in 120 



icit n = 6 calculation of 



21 



221 ] (see also 



111 ] has now demonstrated, following 



411]). that these relations are not true for all 



n, it is believed that they should hold for four and five particles at any number of loops 
since for these cases the amplitudes are determined up to a few constants by dual conformal 
invariance {9], Q] . 
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The precise form of the BDS ansatz at three loops, in dimensional regularization to 
D = 4 - 2e, is 

Mf( e ) = -i(M«( e )) 3 + M«( f )Mf( e ) + /( 3 )(e)M«(3e) + + 0(e) (4.1) 

where 

/(3)(e) = w + (5C(2)C(3) + 6C(5)) e + Ge2 ' c(3) = 6 (4 - 2) 

in terms of two previously undetermined numerical constants a and b. BDS verified by 
explicit calculation that the 3- loop 4-particle amplitude satisfies (14. ip . but the structure of 
the equation for n = 4 is insensitive to the values of a and b as long as they obey the linear 
relation 

341 

2a - 96 = -^C(6) + 17C(3) 2 . (4.3) 

Here we will use our 3-loop 5-particle amplitude to extract a second linear equation from (14. ip 
which will finally fix the constants a and 6. 

The calculation of a and b benefits from two simplifcations. The first is that we may 
restrict our attention to the parity-even part of (14.11) . If we write each amplitude as a sum 
of its parity-even and parity-odd parts, = Me-+ + Mg_ , then taking the parity-even 

part of ( 14. ip for n = 5 gives 

Mff(e) = ~(M«(e)) 3 + M«(^ 

+M«(e) (M£)( e ) _ M«(e)M«(e)) + 0(e). (4.4) 

In Q it was shown that M£ } (e) - M 5 ( + ) (e)M 5 a) (e) = O(e). Since M^(e) itself is also 0(e), 
we see that the entire last line of ( 14. 4p can be replaced simply by +0(e). A consequence of 
the result of is therefore that the parity-even part of the three-loop BDS ansatz can be 
obtained by making the naive replacement Mg -> in dO) . 

The second simplication is to make use of the notion of obstructions introduced in 



and exploited in the four- loop calculations 25|, |38|. We refer the reader to 25j for all of 



24) 



the necessary details, including a detailed algorithm for calculating obstructions. Here we 
simply remind the reader that for an amplitude A(x, e) depending on a single kinematic 
variable x, the obstruction P(e) is defined to be the coefficient of the simple pole at y = 
in the inverse Mellin transform transform, so that 



A(x, e) = / dy x 



Pie) 

(higher order singularities) H h (regular at y — 0) 

y 



14 



(4.5) 



As explained in 25J it is important to understand this relation as holding order by order in e, 
rather than at finite e. The prime advantage of dealing with P(e) rather than the full A(x, e) 
is that it is much simpler to compute. Furthermore it is important that obstructions satisfy 
a product algebra — the obstruction in any product of amplitudes is equal to the product of 
the individual obstructions. 

This concept generalizes straightforwardly to integrals depending on more than one kine- 
matic variable. In the case at hand we have 5-particle integrals depending on five inde- 
pendent variables s^+i, and we can extract the obstruction P(e) by applying the above 
procedure five times in succession. Equivalently, we define -P(e) to be the coefficient of the 
l/(l/il/2l/3l/4l/5) pole in the 5-fold inverse Mellin transform of A(s 12 , s 2 3, s 34 , s 45 , s 5 i). By ap- 



plying the algorithm outlined in 



251 ] it is straightforward to find that the obstructions in 



the one- and two-loop five particle amplitudes are given by 

a) _5 1 5tt 2 17903) 97vr 4 2 _ ^ 51tt 2 C(3) _ 137C(5)\ , 

^ [€) ~ 2e 2 + 8 + 24 e+ 1440 C I 32 8 ]< 



/ 763C(3) 2 _ 23tt 6 \ 
V 72 3780 J 



e 4 + 0(e 5 ) 



,(2), 



, , , 25 1 35tt 2 1 865C 3 1 97tt 4 

py j (e) = r = — + 21.494969e 

5 W 8 e 4 24 e 2 48 e 1152 

64.357473e 2 + C(e 3 ). (4.6) 



For simplicity we have restricted P^ here to the parity-even parts of the amplitudes. Note 
that these expressions satisfy the two-loop ABDK relation 

P?\z) = l(P?\e)) 2 + ("C(2) - C(3)e - C(4)e 2 )P 5 (1) (2e) - ^ + 0(e) (4.7) 
as expected. 

At three loops, we have found that there are nine independent integrals which contribute 
to the parity-even part of the 5-particle amplitude. The obstructions for each of these types 
of integrals, through O(e ), are 

= 20 1 20^ 1 43«3) 1 73^1 _ 1 

9 e 6 27 e 4 2 e 3 432 e 2 e 
70 1 45vr 2 l 1495C(3)1 76vr 4 1 , 1CQnnco . no l 



Pb = ^ - -=--r ~ + 1589.962798- + 2824.770745, 

3 e° 2 e 4 be 6 135 e 2 e 

20 1 25vr 2 1 557C(3) 1 17137vr 4 1 1 
= T 7 -ST 7 ~ ""ST ? + TaST? + 221-894995- + 1030.164974, 
_ f 1 _ 3^ 1 _ 64^ 1 + 1 + ^ m687 i _ 4m 6578ga 

1 155vr 2 1 563C(3) 1 487vr 4 1 1 

P G 2a = 20- - qcc 2 + 1294.520402- + 2938.6610 ± 0.0036, 

e b 9 e 4 4 e J 288 e 

15 



P K = ™L _ _ HIZCW 1 + 11^1 + 178 . 487460 I _ 2387.290195, 

9 e 6 18 e 4 36 e 3 4320 e 2 e 
35 1 85vr 2 1 1411C(3) 1 1195vr 4 1 1 

P u = j — -? -- 673.319831- - 2845.889639, 

3 e 6 12 e 4 12 e 3 432 e 2 e 

1 A^-w 2 1 1 Q837T 4 1 1 

P M = 15-^ - ——-j - 136C(3)4 + ^7^4 + 625.875398- + 437.509754 



e 



36 e 4 e 3 1440 e 2 e 



P R = _80C(3) 1 + IOTttI _ 3g5 562g04 l + g23 415ig6 (4g) 
3 e 3 108 e 2 e 

Each expression displays the result obtained after summing over all 10 permutations of the 
corresponding integral (including in each case the appropriate dual conformal numerator). 
The estimated error in the numerical results is much smaller than the precision indicated in 
all cases except for the last term in Pc2a, which is the overwhelmingly dominant source of 
numerical error. 

Using the parity-even parts of the coefficients obtained in the previous section, and includ- 
ing the necessary factors of 1/2 to avoid overcounting those integrals with flip symmetries, 
we find the total three-loop obstruction 

-P 5 (3) = ( p Ai - ^Pb + Pd + ^Pd 2 + Pc2a + Pk + -P^ + f/Vj - 2 P R^j 

125 1 325vr 2 1 4175C(3) 1 499vr 4 1 
~ ~l8~^ + 192 ? + 192 1? + 10368^ 

-40.764885- + 207.1613 ± 0.0002 + 0(e) (4.9) 
e 

Using the results (14.91) and (14.61) . we find that the BDS relation (14.11) is satisfied provided 
that a and b satisfy the linear relation 

5a - 186 = 105.482 ± 0.004. (4.10) 

Together with (14.31) this implies the solution 

a = 85.263 ± 0.004, b = 17.8241 ± 0.0009. (4.11) 

Based on the transcendentality hypothesis, it is expected that each of these numbers should 
be a linear combination of £(6) and C(3) 2 with rational coefficients. However given the limited 
numerical accuracy of our calculation it seems prudent to avoid speculating on possible exact 
values for a and b at this time. 
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V. SUMMARY 



In this paper we have used the leading singularity method to obtain an ansatz for the 
four- dimensional cut-constructible part of the 3-loop 5-particle amplitude in M = 4 SYM 
theory. This means that we have determined the coefficients of the integrals shown in fig. [2] 
by comparing residues of the ansatz to those of the amplitude on various leading singular- 
ities. Although it has not yet been proven that determination of only leading singularities 
completely determines an amplitude (in principle one might have to add additional integrals 
that vanish on all leading singularities but that have subleading singularities), the method so 
far has been found to give the complete answer in all cases where comparison with alternate 
methods was possible. 

Dimensionally regulated amplitudes occasionally contain so-called '//-terms' which are 
defined as terms in the integrand which vanish in D = 4 but not in D = 4 — 2e (note that 
this statement is, in general, completely unrelated to whether or not these terms vanish 
in D = 4 — 2e after integration; indeed //-terms can easily be IR divergent). Since the 
leading singularity method itself works with strictly four-dimensional loop momenta, it is 
insensitive to possible // terms, although it seems that in principle they could be determined 
by considering leading singularities in integer dimensions other than 4. However, in all cases 
that have been studied so far it has been observed that //-terms separately cancel out of 
the BDS relation, leaving unaffected. We can therefore hope that even if the 3-loop 
5-particle amplitude contains such terms which we have missed, they would not contribute 
to the constants a and b computed in section IV. 

Finally we emphasize that since our goal in section IV was to streamline the calculation of 
a and b as much as possible, we have only evaluated the obstructions, not the full amplitude. 
Consequently we have not checked (even numerically) that the quantity appearing 
in (I2.3P is a numerical constant; in principle it could depend on the kinematic variables 
s^j+i. The method of obstructions is efficient for quickly extracting the 'constant part' of 
(defined as the coefficient of 1/y in the inverse Mellin transform) but is insensitive to 
any other potential terms in that depend on the s^j+i. It remains an interesting open 
problem to verify that there are no such terms. 
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